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We present sub-millimeter statistical detections of galaxies discovered in the 
5' x 5' Spitzer Early Release Observations (to ~ 4 — 15/xJy 5<r at 3.6 — 8/xm, 
00 ' 170/iJy at 24/im) through a stacking analysis of our reanalysed SCUBA 8mJy 

survey maps, and a Spitzer identification of a new sub-millimeter point source 
in the 8mJy survey region. For sources detected at 5.8/xm or 8/im (154 and 



O 



111 sources respectively), we detect positive skews in the sub-millimeter flux 
distributions at 99.2 — 99.8% confidence using Kolmogorov-Smirnov tests, at 
both 850/im and 450/mi. We also marginally detect the Spitzer 24/im galaxies at 
850/im at 97% confidence, and place limits on the mean sub-millimeter fluxes of 
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the 3.6/zm and 4.5/zm sources. Integrating the sub-millimeter fluxes of the Spitzer 
populations, we find the 5.8/im galaxies contribute 0.12 ± 0.05 nW m -2 sr~ 4 to 
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the 850/xm background, and 2.4 ± 0.7 nW m~ 2 sr _1 to the 450/im background; 
similar contributions are made by the 8/xm-selected sample. We infer that the 
populations dominating the 5.8/im and 8/xm extragalactic background light also 
contribute around a quarter of the 850yum background and the majority of the 
450/im background. 

Subject headings: cosmology: observations - - galaxies: evolution - - galaxies: 
formation 



1. Introduction 

The early pioneering sub-millimeter surveys in lensing clusters (Smail, Ivison & Blain 
1997) and in blank fields (Hughes et al. 1998, Barger et al. 1998) demonstrated the feasibil- 
ity of deep extragalactic surveys exploiting the favorable K-corrections in the sub-millimeter. 
Deep 850/zm imaging has now resolved around half of the 850/iin extragalactic background 
(Hughes et al. 1998, Blain et al. 1999, Cowie, Barger & Kneib 2002). These galaxies are of- 
ten called SCUBA galaxies after the instrument with which they were first detected (Holland 
et al. 1999). The 14" 850/zm SCUBA beam makes identifications at other wavelengths dif- 
ficult; nevertheless, ~ 50% of 850/mi sources are identifiable in /zJy 1.4GHz imaging (Ivison 
et al. 2002). These radio identifications have led to optical identifications, morphologies and 
ultimately spectroscopic redshifts in multiwavelength follow-up campaigns (e.g. Chapman 
et al. 2003). Furthermore, the non-detection of SCUBA galaxies in hard X-ray imaging 
(e.g. Alexander et al. 2003) suggests that the bulk of the population has far-infrared lumi- 
nosities dominated by star formation. The morphologies, redshifts, clustering and molecular 
gas contents are so far consistent with at least some of the SCUBA population being the 
progenitors of giant elliptical galaxies (e.g. Dunlop 2002), though other alternatives are still 
viable (Efstathiou & Rowan-Robinson 2003) and the SCUBA population is heterogeneous 
(e.g. Ivison et al. 1998, 2000). Finally, the K-correction effects in the sub- millimeter make 
the sub-millimeter extragalactic background sensitive to contributions from the far-infrared 
luminous energy densities at all redshifts z < 10. The populations which contribute to the 
850/im extragalactic background are necessarily also significant contributors to the cosmic 
history of dust-shrouded star formation. 

Following the IRAS mid-infrared surveys of the local Universe (e.g. Rush, Malkan & 
Spinoglio 1993), the mid-infrared was first made accessible to deep extragalactic surveys by 
the Infrared Space Observatory (ISO, Kessler et al. 1996) which conducted a suite of surveys 
with a variety of depths and areal coverages (e.g. Genzel & Cesarsky 2000 and refs. therein, 
Rowan- Robinson et al. 2004). The rapid upturn in the 15/im extragalactic source counts 
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clearly demonstrated the existence of a strongly evolving population of obscured starbursts 
and active galaxies (e.g. Serjeant et al. 2000, Franceschini et al. 2001, Gruppioni et al. 2002). 
It has also been argued that the populations dominating the 15/zm extragalactic background 
light, which are resolved by ISO, are also largely the same populations which dominate 
the unresolved 140/zm background (Elbaz et al. 2002). If correct, this is a significant 
breakthrough in determining the populations which dominated the far-infrared luminous 
energy density throughout the relatively recent history of the Universe (e.g. z < 1). 

Finding the population that supplies the luminous energy density at z > 1 requires 
understanding the sub- millimeter background light. However, it has been difficult to find 
sub-millimeter source counterparts in the mid-infrared. Very few sub-millimeter-selected 
sources have been detected by ISO in the mid-infrared (e.g. Eales et al. 2000, Webb et al. 
2003b, Sato et al. 2002). The reverse procedure of looking for mid-infrared sources in the 
sub-millimeter via stacking analyses have not fared much better. Serjeant et al. (2003a) 
found no excess 850/im flux at the locations of 15yum sources in the HDF North. Lyman 
break galaxies, in contrast, are detectable statistically (e.g. Peacock et al. 2000, Webb et 
al. 2003a). If SCUBA galaxies are extreme star-forming galaxies in the most massive high- 
redshift halos, then their anomalously faint K-band identifications imply heavy obscuration 
in the observed-frame near- infrared (Serjeant et al. 2003b), suggesting that SCUBA galaxies 
may be detectable in /zJy-level mid-infrared imaging. 

The Spitzer Space Telescope (hereafter Spitzer, Werner et al. 2004) is an enormous 
advance over ISO in terms of mapping speed, sensitivity, and wavelength coverage. It may 
now be possible to resolve the bulk of the extragalactic background light at 3.6 — 8/xm in 
exposures of the order 15 minutes, an equivalent depth achieved at 6.7/zm in 23 hours with 
ISO by Sato et al. 2003. In this paper we present statistical sub-millimeter detections of 
galaxies selected in a new Spitzer survey. The Spitzer identifications of previously-published 
sub- millimeter sources are discussed by Egami et al. (2004). Identifications of MAMBO 
sources are discussed by Ivison et al. (2004). 



2. Observations 

The Spitzer Early Release Observations survey is one of the first extragalactic surveys 
conducted by Spitzer. Besides the pragmatic goal of characterising the survey capabilities of 
the facility, the survey has the science goals of making the first constraints on the populations 
which dominate the extragalactic backgrounds in the shorter wavelength Spitzer bands, and 
the links between these galaxies and other known populations. Accordingly, the survey field 
was selected to lie in the Lockman Hole, an area with abundant multi-wavelength survey 
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coverage, and in particular with 7 galaxies from the 850/im 8mJy survey > 3.5a catalog 
(Scott et al. 2002, Fox et al. 2002, Ivison et al. 2002; see below). 

The Spitzer imaging is described by Egami et al. (2004) and Huang et al. (2004). In 
summary, IRAC (Fazio et al. 2004) imaged a ~ 5' x 5' field for 750s per sky pixel at all four 
bands, resulting in median la depths of 0.77/iJy, 0.78/iJy, 2.75/iJy and 1.67/iJy at 3.6/im, 
4.5/im, 5.8/im and 8/im respectively, and sources were extracted to > 4.7o\ MIPS (Rieke et 
al. 2004) observed the field in photometry mode at 24/im for 250s per sky pixel, resulting in 
a typical la depth of 30/iJy. Source confusion and blending make the MIPS completeness 
and reliability more problematic than for SCUBA or IRAC. MIPS source catalogs were 
extracted to > 170/iJy (80% complete, > 90% reliable) and to > 120/iJy (50% complete, 50% 
reliable). We conservatively assume the catalogs are all extragalactic, as expected at these 
flux densities; the deletion of any (hypothetical) contaminant population of galactic stars 
with negligible sub-millimeter fluxes would improve the confidence levels of our statistical 
detections. 

The SCUBA 8mJy survey data (which covers all our Spitzer field) is described in detail 
by Scott et al. (2002). We have reanalysed this data, modelling the intra-night SCUBA gain 
variations, improving the extinction corrections using improved fits to the 225GHz skydips 
monitored at the Caltech Sub-millimeter Observatory, and removing the cross-talking 450/im 
bolometers A7 and A16 (see Mortier et al. 2004). Additional sources are detectable in the 
zero-sum chopped/nodded sub-millimeter maps. After point source detection using noise- 
weighted PSF convolution (Serjeant et al. 2003a) including the effect of negative chop/nod 
positions, the median la depth (total flux) in the Spitzer Early Release Observations field 
is 1.6mJy at 850/im, and 14.4mJy at 450/im. 



3. Results 

Our reanalysis of the SCUBA 8mJy survey uncovered new candidate sub-millimeter 
sources (Mortier et al. 2004). Figure 1 shows the Spitzer identifications of the most significant 
of these, at a = 10 h 51 m 53.94 s , 5 = +57° 25' 05.5" (J2000), detected with a 3.7a 850/im flux 
of 4.3mJy. At 450/im the source is not detected, though the 450/im flux at its position is 
17.9mJy (1.6<j); there is also a hint of positive flux at 1.8c significance at 1250/im (Greve et 
al. 2004) at the position of this source. There are hints that the source is extended at 850/im, 
suggesting a blend of more than one source, and indeed there are two candidate identifications 
in the Spitzer imaging. The probabilities of a random association (P = 1 — exp(— n(> s)-Kr 2 ) 
where n(> s) is the surface density of catalogued objects brighter than the identification 
at distance r) are 0.033 and 0.155 for the brighter and fainter 4.5/xm flux respectively. 
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Interestingly, the low-significance contours of 450/im emission are coincident with the brighter 
of the two candidates. The Spitzer fluxes of this source are 16.6, 26.7, 14.9, 18.9/iJy and 
216/iJy at 3.6/iin, 4.5/xm, 5.8/im, 8/xm and 24/im respectively. The fluxes of the second 
candidate identification are 6.2/iJy and 4.6/iJy at 3.6//m, and 4.5/im, but the source is 
not detected at longer wavelength Spitzer bands. While deblending the fluxes from the 
two identifications is beyond the scope of this paper, it is interesting to note that either 
identification would have a sub-millimeter: Spitzer flux ratio which is redder than in the 
population as a whole, derived below. 

The > 3.5a sub-millimeter sources are already known to have Spitzer identifications 
(Egami et al. 2004), so we construct a mask to remove the sources at these positions from 
the stacking analysis. At 850/im we masked a 9.9" radius (V2x the SCUBA FWHM, i.e. the 
FWHM of point sources in our PSF-convolved maps) around the > 3.5a 850/im point sources 
from Scott et al. (2002), as well as the new point sources from our reanalysis (Mortier et al. 
2004). At 450/im the masking radius was 5.2", which we also applied to new 450/im sources 
from our reanalysis. We also masked regions with high noise levels, arbitrarily selected as 
> 5mJy at 850/im, or > lOOmJy at 450/im, and also masked all regions in the sub-millimeter 
maps without Spitzer data. The unmasked area at 850/im (450/im) is 22.1 (22.7) arcmin 2 
for the IRAC catalogues. For MIPS, the corresponding area is 27.6 (28.5) arcmin 2 . 

The sub-millimeter maps give the best-fit point source flux at each location (equation 
A4 of Serjeant et al. 2003a), so we measure the values of the sub-millimeter images at the 
Spitzer galaxy positions. Figure 2 shows the sub-millimeter signal-to-noise ratios (equation 
A6 of Serjeant et al. 2003a) at the positions of the Spitzer Early Release Observations 
source catalogs at 5.8/im and 8/iin, which lie in the unmasked regions. These figures also 
show histograms for the whole of the unmasked maps. Note the clear positive skews in the 
sub-millimeter fluxes at the Spitzer source positions, relative to the maps as a whole. We 
used Kolmogorov-Smirnov tests to determine the confidence level of these relative positive 
skews, listed in table 1. There are significant detections at 5.8/im and 8/iin, and marginal 
detections at 24/im. No significant skew was found at 3.6/im or 4.5/im. 

We also calculated the mean fluxes in the sub-millimeter maps at the positions of the 
Spitzer galaxies. These are almost all positive and are listed in table 1. This statistic is less 
efficient than the Kolmogorov-Smirnov test (which does not just use the first moment), so 
the uncertainties on the mean fluxes are larger than suggested by the Kolmogorov-Smirnov 
confidence levels. The Spitzer galaxies are individually undetected in the sub-millimeter, 
so one must be careful in interpreting the stacking analysis; for example, we do not use 
the noise-weighted mean fluxes because they can give false positives, as noted by Serjeant 
et al. 2003a. Also, many fluctuations in the sub-millimeter maps are blends (e.g. Scott 
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et al. 2002), so is there a risk of overestimating the flux of any given Spitzer galaxy by 
also counting its neighbors? The answer is no, as Peacock et al. (2000) showed: for a 
Poissonian sampling of the sub-millimeter map, the expected total sub-millimeter flux from 
all neighbors equals the mean flux of the map, which is exactly zero for the chopped, nodded 
SCUBA maps. (The chopped/nodded PSF is also zero-sum, so the PSF-convolved images 
are still zero-sum.) Blending slightly degrades the uncertainty on the mean, but it does not 
affect the mean fluxes. Another demonstration of this is to note that the total flux of any 
given source is exactly zero in these chopped, nodded maps. Even confusing sources within 
the same SCUBA beam will have no net effect, provided that there are also similar sources 
uniformly distributed over the rest of the map. 

The histograms of the whole map in figure 2 effectively act as the control sample. Nev- 
ertheless, as a test of the stacking analysis methodology, we performed the same tests on 
simulated Spitzer catalogs of the same size as the observed catalogs, with a Poissonian dis- 
tribution in the unmasked regions. No positive skews in the sub-millimeter flux distributions 
were detected relative to the maps as a whole. These simulations also verified that the mean 
sub-millimeter flux around randomly chosen positions is consistent with zero. 



4. Discussion 

Although our Spitzer 5.8/im and 8/im catalogues are small (154 and 111 respectively), 
we can use our statistical detections to obtain constraints on the contribution to the sub- 
millimeter extragalactic background light from the Spitzer populations, by multiplying the 
mean sub-millimeter flux of a single Spitzer galaxy (table 1) with the observed surface density 
of Spitzer galaxies. At 850/im, the integrated contribution from the 5.8/im population is 
0.12 ± 0.05 nW m~ 2 sr _1 , and from the (overlapping) 8/im population 0.10 ± 0.04 nW m~ 2 
sr _1 . This is already comparable to the total 850/im extragalactic background reported by 
Lagache et al. (2000) of ~ 0.5 nW m~ 2 sr _1 . At 450/im the Spitzer contributions can, within 
the uncertainties, account for all the observed ~ 3 nW m -2 sr _1 extragalactic background: 
2.4 ± 0.7 nW m~ 2 sr _1 from the 5.8/iin-selected sample, and 2.2 ± 0.6 nW m~ 2 sr _1 from the 
8/im galaxies. Our sub-millimeter faint source density is consistent with those in the deepest 
850/im maps (Cowie et al. 2002). 

The sub-millimeter:>S'piteerflux ratios in sub-millimeter point sources (Egami et al. 2004, 
see above) are much higher than average in the population (table 1). Despite the common- 
ality between the SCUBA and Spitzer galaxy populations, sub-millimeter galaxies are not at 
all representative of the Spitzer population. The sub-millimeter background independently 
suggests the same result: if the 24/im population had sub- millimeter -.Spitzer flux ratios as 
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high as SCUBA galaxies (e.g. median S^o^m/S^m = 34.4 for the Egami et al. (2004) 
sample), the > 170/zJy 24/xm sources alone would over-predict the 850/iin background by 
~ 60%. Our weak sub- millimeter detection of the 24/zm population is less surprising in 
that context. One possible interpretation is that the SCUBA population has fewer (or more 
heavily obscured) active nuclei than the Spitzer population as a whole. 

Is our stacking analysis signal due to dusty galaxies, or to high- z galaxies, or both? 
Our samples are too small to distinguish these possibilities, but stacking analyses of subsets 
show that either option is credible. Almost all the 5.8/xm galaxies are detected at 3.6/im, 
and as with the 5.8/xm population as whole (table 1) this subset (153 galaxies) is detected 
statistically. The 299 remaining 3.6/xm galaxies are not detected in the stacking analyses 
(e.g., only 37% confidence in the 450yum signahnoise maps), so the redder galaxies are more 
likely to be sub-millimeter emitters. The lack of a stacking signal at 3.6/xm and 4.5/iin is at 
least in part due to dilution of the signal from the high density of bluer galaxies detected 
at these wavelengths. The red colors could be due to heavy obscuration, but the colors 
also become redder quickly with increasing redshift. Simpson & Eisenhardt (1999) show 
that i 7 4. 5/im /F 3 . 6/im > 1.32 or F 8fim /F 5 .8^ m > 1.32 selects redshifts z > 1.5. There are 47 
galaxies satisfying either criteria, with detections in the relevant bands, and these galaxies 
are detected in the stacking analysis (e.g. 99.1% confidence at 450/xm); the corresponding 400 
foreground z ;$ 1.5 galaxies are not detected (e.g. only 71% confidence at 450yum). Further 
work on larger samples is needed to discriminate between high-redshift and high-obscuration 
sub-populations. 

The sub-millimeter backgrounds are sensitive to contributions from the far-infrared lu- 
minous energy densities throughout the redshift range of favorable K-corrections. It follows 
that the Spitzer 5.8/xm and 8/im-selected galaxies are necessarily significant contributors to 
the comoving volume averaged star formation density. Many of the 5.8/xm and 8/xm galaxies 
are challenging targets for 8-10m-class optical/near-infrared spectroscopy, but redshift sur- 
veys of these populations, together with 2-D and 3-D clustering, would discriminate between 
competing semi-analytic descriptions of galaxy evolution, and determine whether (or which) 
SCUBA/ Spitzer galaxies are the sites of the assembly of giant ellipticals. Finally, it is worth 
stressing that the conclusions in this paper could only have been reached with co-ordinated 
mult i- wavelength surveys. 
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Fig. 1. — Spitzer identification of the sub-miilimeter source in section 3. All postage stamps are 
30" on a side, with North up and East to the left. From left to right, the greyscale images show the 
Subaru R-band image; 4.5//m; 5.8/mi and VLA 1.4GHz A-array. The contours, from left to right, 
are 850/um, 450/um, 24/um and 850^111. In all cases the contours are from la in steps of 0.5o\ 
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Fig. 2. — Histograms of the sub-millimeter signal-to- noise ratios (see text) at the positions of the 
Spitzer 5.8/um and 8/xm galaxies (hatched histograms, left-hand ordinates), compared to the his- 
tograms for the unmasked regions of the sub-millimeter maps as a whole (unhatched histograms, 
right-hand ordinates). Note the clear positive skew in all four hatched histograms, relative to the 
unhatched histograms. The diagrams are annotated with the probabilities that identical distribu- 
tions would differ by as much as is observed, according to Kolmogorov-Smirnov tests (see also table 
1). 
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Sub-millimeter stacking analysis of Spitzer source catalogs 
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lable 1 '. Stacking analysis of the Spitzer galaxy populations in the sub-millimeter maps. The 24^im > 120^tJy and > 170/iJy 
catalogs are listed separately ( a and " respectively). The N svc values are the numbers of Spitzer sources used in each stacking 
analysis. The Kolmogorov-Smirnov test results quoted are confidence levels, i.e. a high number indicates a net statistical 
detection of Spitzer galaxies. Results are listed for the signal-to-noise histograms (plotted in figure 2) as well as the flux 
distributions. At 5.8^im and 8fim the galaxies are detected, marginally so at 24^tm, but not at 3.6fj,m or 4.5/im. Also given are 
mean sub-millimeter fluxes, and mean flux ratios, at the positions of the Spitzer sources. Errors quoted are a/yN, where N 
is the number of measurements and cr 2 the variance of the measurements. Known sub-millimeter sources have been excluded 
from this analysis, as discussed in the text. 



